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1 .  Introduction 

Hln  a  performance  comparison  of  reentry  vehicles  there  are 
a  number  of  important  interrelated  criteria  determining  the  ac¬ 
ceptability  of  the  vehicle.  These  include  the  peak  deceleration, 
the  limits  imposed  by  atmospheric  heating,  guidance  control  and 
accuracy,  and  time  of  descent.  A  general  study  has  been  under¬ 
taken  which  places  emphasis  on  the  first  two  considerations  and 
is  based  on  three  vehicle  configurations  appropriate  to  both 
manned  and  unmanned  entry  into  the  earth' s  atmosphere  from  high 
supercircular  speeds.  The  first  vehicle  type,  designated  "sphere" 
is  a  high  drag,  nonlifting  body.  The  second,  referred  to  below 
as  the  "body"  type,  is  of  moderate  drag  and  possesses  a  trimming 
lift  capability;  while  the  third,  the  "wing",  simulates  low-drag, 
high  lift  vehicles.  For  each  configurational  type  heat  transfer 
behavior  is  to  be  assessed  on  the  basis  of  three  techniques  of 
heat  dissipation,  viz.,  by  radiation  alone,  by  ablation  or 
internal  cooling,  and  by  the  utilization  of  the  vehicle  as  a 
heat  sink. 

In  this  paper  are  presented  the  results  of  the  heat  sink 
studies  in  the  form  of  thermal-time  response  and  specifically 
in  terms  of  temperature-time  histories  and  peak  temperatures  as 
a  function  of  vehicle  mass,  shape  and  aerodynamic  parameters, 
entry  orbit.  The  thermal  response  is  correlated  to  peak 
deceleration. (  ^ 

The  nose  of  each  vehicle  can  be  thought  of  as  spherical 
and  is  defined  by  the  nose  radius  R  and  the  wall  thickness  D. 

t 

If  R  and  D  are  within  certain  bounds,  the  differential  equation 
governing  the  temperature  of  the  nose  is  essentially  one-dimen¬ 
sional.  The  heat  balance  is  dependent  on  the  velocity  and 
altitude  of  the  vehicle  at  all  times.  These  are  known  from 
earlier  trajectory  computations.  For  a  representative  range 

of  single  and  multipass  entry  orbits  there  are  presented  peak 

• 

temperature  contours  in  the  R  -D  plane  which  permit  a  choice 
of  R  ,  D  combinations  which  for  a  given  peak  deceleration  place 
limits  on  temperature. 


2 


In  determining  material  constants  appearing  in  the 
differential  equation  the  heat  shield  is  assuwned  to  be  manu¬ 
factured  of  a  so-called  high  performance  material  such  as 
beryllium  oxide  (T  =5125°R). 

m  a  X 

2.  Notation 

a'  =  25,688 

B*  =  1349.2 

Cq  =  -7.348}  Cj  =  .10842}  C2  =  -3.1476.10~5 

C3  =  56.25}  C4  =  .0175 

D  =  0.01  -  0.05  feet}  thickness  parameter,  variable 

e'  =  0.21625. 10~12  .  T*  =  .055303  (BTU/ft2sec) 

e  eccentricity  of  ex-atmospheric  ellipse 

h  time  step  used  in  temperature  calculations.  (sec) 

I  I 

h  time  step  used  in  trajectory  calculations.  (sec)(h=h  /N) 

k  =  0.11856  •  109  ft.3/2  sec”1  . 

N  =  h'/h 

« 

R  =2  ft,  -  12  ft.}  variable  vehicle  parameter. 

R  radial  distance  from  the  earth's  center}  non-dimensional ized 

with  respect  to  r  . 

O 

Rq  =  .96957  =  radius  of  spherical  earth 

r  =  2.1556  •  107  feet}  radius  of  the  initial  condition 

3 

(200  km. ) 

2  2 

S  =2  cal/min.  cm.  =  .1229  BTU/sec.ft.  }  solar  constant 

I  I 

t  =  t  (8)  time  in  seconds  in  an  elliptical  orbit,  measured 

from  the  line  of  apses. 

•I  M 

t  =  t  (T)  time  in  seconds  (see  "Multiple  Pass  Trajectories"). 

T  =  T/T  }  non-dimensionalized  wall  temperature  of  vehicle, 

w  s 

T  absolute  wall  temperature  of  vehicle. 

Ts  =  711.13°R  =  (aS/E  J1/^  equilibrium  temperature  of  the 

vehicle  outside  the  atmosphere. 
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Y  non-dimensional  velocity  =  v/vq 

VQ  initial  velocity  V  at  200  km.  altitude  for  which  R=l. 
v  absolute  velocity  (ft. /sec.). 

vq  25,546  ft. /sec.}  circular  orbital  velocity  at  200  km. 
a  It i tude  . 

v  =  26,000  ft. /sec.}  normalizing  factor  for  velocity  used 

s 

in  temperature  analysis, 
y  altitude  above  sea-level  (feet). 

ye  =  0.65595. 106  ft.}  initial  altitude  (200  km.), 
a  =  0.45}  vehicle  coefficient  of  absorptivity, 

a  angle  of  incidence 

p  =  1/22,000  =  .45455  10"4  ft.-1  scale  height  of  the  atmos¬ 

phere. 

T  =  period  of  elliptical  orbit 

u  angular  coordinate}  true  anomaly  in  Keplerian  ellipse. 

Y  flight  path  angle,  positive  upwards 

Y  initial  value  of  y, 

'  o 


3.  The  Heat  Transfer  Equation 

The  one-dimensional  differential  equation  governing  the 
non-dimensional  wall  temperature  Tw  of  the  vehicle  isi 


dT 


w 


dt 


=  t 


(R  ) 


172 


(— >13/4[l-0.3  T„/hse]e-^/2-E' 


1 

C(TW)D 


This  states  that  the  heat  transfer  is  a  balance  of  forced 
convection  heating  (aerodynamic  effects),  thermal  radiation 
(according  to  Planck's  law)  and  solar  irradiation.  The  source 
of  the  first  term  in  this  equation  is  Ref.  (l),  and  it  is  seen 
to  include  the  influence  of  surface  temperature  on  the  convective 
heating  rate.  A  value  of  the  solar  constant  S  was  found  in 
Ref.  (2).  Other  numerical  data  are  listed  in  the  Notation. 

The  coefficients  of  the  differential  equation  depend  on 
the  time  history  of  velocity  v  and  altitude  y}  and  can  be  computed 
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from  pre-cal cul a  ted  trajectories  given  in  the  form  of  V,  R 
( nondimensional  velocity,  radial  distance)  pairs  at  discrete 
time  intervals. 


4.  The  Numerical  Integration 

The  temperature  equation  can  be  integrated  numerically 
by  a  Runge-Kutta  technique  (Ref.  3).  In  the  procedure,  an 
equation  of  the  type  dT/dt  =  f(t,  T)  is  integrated  using  the 
well-known  formulasj 

T ( t+h)  =  T ( t)  l/6(k1+2kp+2k3+k4)  , 

kx  =  hf(t,T);  k2  =  hf  ( t+h/ 2 ,  T+^/2)  , 

k3  =  hf ( t+h/ 2 ,  T+k2/2){  k4  =  hf(t+h,  T+k3). 

The  time  step  h  must  be  chosen  sufficiently  small.  The  per¬ 
missible  value  is  determined  by  successively  decreasing  h  until 
results  become  insensitive  to  variations  in  h. 

In  the  present  problem,  the  velocities  V  and  altitudes 
y  (or  R)  which  determine  f(t,T),  are  tabulated  at  time  intervals 
h  .  The  time  interval  used  in  the  trajectory  calculations,  h  , 
usually  varies  considerably  during  any  given  trajectory  (Ref.  4). 
In  the  temperature  calculations,  the  time  step  h  was  determined 
from  the  given  time  step  h  by  the  relation  h=h  /N,  N  being  an 
integer.  Thus  for  N=1  ,  the  temperature  is  calculated  at  the 
same  time  intervals  as  the  given  velocity  and  altitude.  The 
corresponding  velocity  and  altitude  are  interpolated  linearly 
between  any  two  values.  It  was  found  by  comparative  calculations 
that  parabolic  or  higher  order  interpolation  between  the  table 
values  was  unnecessary. 

The  Runge-Kutta  process  was  programmed  for  a  Bendix  G15 
Computer,  so  that  there  are  printed  out  six  quantities  for  each 
integration  step  in  the  following  orderj  (t)  time  in  seconds; 

(T  )  wall  temperature  nondimensional  i zed  with  respect  to  T&  = 
711.13  R;  (T)  temperature  in  °R;  (V)  given  nondimensional  i zed 
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velocityj  (R)  given  nondimensional ized  distance  from  the  earth's 
center*  (dT/dt)  =  f(t,  T)  rate  of  change  of  nondimensionalized 
temperature.  The  time  required  to  compute  one  temperature  is 
approximately  50  seconds. 

Comparative  results  indicate  that  the  time  step,  h,  used 
in  the  temperature  calculations  can  be  chosen  equal  to  the  given 
time  step,  h  ,  of  the  trajectory  calculations  which  is  related 
to  the  acceleration  magnitude*  i.e.,  N=1  is  acceptable.  The 
initial  temperature  of  the  vehicle,  which  must  be  prescribed, 
can  be  varied  by  as  much  as  10#  without  noticeably  influencing 
the  ensuing  temperature  history.  The  initial  temperature  (i.e., 
the  temperature  Tg  =  (2S/E  ) of  a  body  in  free  space  for 
which  thermal  equilibrium  exists  between  solar  irradiation  and 
thermal  reradiation)  for  all  atmospheric  passes  was  therefore 
fixed  at  Tg  =  711.13°R. 


5.  Initial  Conditions  and  Input  Data 
Initial  Temperature* 

For  all  calculations  the  chosen  initial  temperature  was 
711.13°R,  which  is  based  on  an  absorptivity  a  =  0.45. 

Time  step* 

The  chosen  time  step  was  the  one  used  in  the  trajectory 

(4) 

calculations*  this  varied  continuously  from  50  seconds  prior 
to  actual  entry,  through  5-10  seconds  near  entry, to  a  minimum 
of  2.5  seconds  when  decelerations  are  greatest. 

I 

R  (nose  radius)  -  D  (thickness)i 

f 

Each  trajectory  was  run  for  several  R  and  D  combinations 

f 

in  the  range  2  <  R  <6  and  .01  <  D  <  .05  ft.  The  combinations 

were  chosen  so  as  to  delineate  a  peak  temperature  contour  in  the 
*  0 

R  -  D  plane  of  close  to  5125  R,  the  acceptable  limit  for  a 
material  such  as  beryllium  oxide. 

"Sphere"  -  trajectories* 

These  trajectories  carry  the  notation  "S".  All  trajec¬ 
tories  begin  outside  the  effective  limit  of  the  earth's  atmos- 


■**■-*'" .  *•»*•«  t 
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phere,  i.e.,  at  an  initial  value  R=1 .  Two  initial  super-circular 
velocities  were  considered*  V=1.3911  and  V=1.3088.  This  group 
includes  nine  single-pass  trajectories,  two  double-pass,  one 
triple-pass  and  one  quadruple  pass  trajectory  defined  by  a 
range  of  entry-angles. 

"Body”  -  trajectories* 

These  trajectories  carry  the  notation  MBM.  These  also 
have  an  initial  R=1 ,  and  an  initial  velocity  V=1.3911.  This 
group  is  composed  of  fifteen  trajectories  including  three  skip 
trajectories  of  which  nine  lead  to  impact  upon  the  first  pass. 
"Wing*'  -  trajectories* 

No  temperature  calculations  have  been  made  so  far. 


6.  Multiple  Pass  Trajectories 

In  order  to  cater  for  multiple-pass  trajectories,  the 
machine  storage  of  the  given  trajectory  data  V,  R,  t  becomes 
excessive.  It  proves  unnecessary  anyway.  Since  the  convec¬ 
tion  terms  in  the  differential  equation  can  be  neglected 
outside  the  atmosphere,  the  equation  can  be  integrated  analytic¬ 
ally.  The  result  is* 

T  +1 

t"(T)  =  (D/4E'Tg){(ao+a2Tg)log  (-!~-)+2(ao-a2Tg)tan"1(T) 

w 

t2+i  « 

+  a, T  log  (-f — )  -  t  (1700)}  , 
l  s  e  T^_1 

w 

where 

(ao,ai,a2)  =  ( Co ,Ci ,C2)  for  500°  <  T  <  1700° 

(a0,a1,a2)  =  (C3,C4,0)  for  1700°  >  T  . 


.. 

.  '  1  ( 

I 

» 

The  time  t  (u)  along  an  elliptical  orbit  of  eccentricity 
e  and  having  a  period  t  is  (Ref.  5) t 

,  1/2  ,  1/2 

t'  =  t  - Ce  5l"  u  +  - - tan(u/2)  ]} 

2it  *l+e  cosu  (1-e  ;  (1+e) 

It  is  Uv t  possible  to  give  an  explicit  relationship  for 

I 

the  temperature  in  terms  of  exit  angle.  A  plot  of  both  t  (u) 

ff 

and  t  (T)  (for  D=.02)  is  given  with  which  it  is  possible  to 
determine  the  temperature  of  the  vehicle  at  entry  into  the 
atmosphere  at  the  end  of  an  elliptical  pass  in  free  space. 

The  information  required  for  this  is  the  temperature  at  exit, 

(i.e.,  at  the  beginning  of  the  elliptical  path)  the  eccentricity 
and  the  period  of  revolution  of  the  ellipse.  See  Figs.  1,  2 
and  3  and  the  Appendix. 

A  study  of  these  graphs  shows  that  for  all  practical 
purposes,  the  satellite  vehicle  cools  off  to  essentially  the 
equilibrium  temperature,  T  =  711.13°R,  during  its  flight 
along  the  ellipse,  irrespective  of  the  temperature  it  had 
at  the  beginning  of  the  ellipse.  However,  for  flight  times 
outside  the  atmosphere  of  less  than  about  10,000  sec.  (for 
D=.02)  the  temperature  at  reentry  can  be  considerably  higher 
than  711.1 3°P . 


7 .  Presentation  of  Results 

Two  forms  of  graphical  presentation  are  given.  The  first 
(Figs.  4-10)  consists  of  wall  temperature-time  histories  and  is 
designed  to  illustrate  the  typical  effects  of  entry  with  and 
without  lift,  of  initial  velocity  (VQ)  and  entry  angle  ( y0 ) > 
and  of  the  loading  ratio  C^S/m.  Temperature  effects  are  related 
to  trajectory  type  and  the  acceleration-time  history  or  the  peak 
deceleration.  Figure  5  relates  the  heat  transfer  (dT/dt)  to 
temperature. 


■4 
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The  second  set  of  curves  (Figs.  11-18)  depicts  the  influ- 

t 

ence  of  vehicle  nose  radius  (R  )  and  skin  thickness  (D)  on  peak 
temperature  for  a  variety  of  initial  velocities,  entry  angles, 
C-S/m  ratios,  and  angles  of  incidence.  Contours  of  constant 

f 

peak  wall  temperature  are  given  in  the  R  -D  plane.  These  were 
obtained  by  linear  interpolation  of  peak  temperatures  at  points 
of  a  mesh,  some  of  which  are  shown  in  addi+ion  to  the  interpolated 
contours. 

Some  essential  properties  of  the  trajectories  used  for 
the  heat  transfer  computations  are  summarized  in  Table  I  for 
the  nonlifting  vehicle  and  in  Table  II  for  the  "body"  with 
trimming  lift  capability.  The  latter  is  defined  by  the  coupled 
lift-drag  laws 


C 


L 


=  Mo 


C 


D 


+  Qa 


2 


where  M  =  QC^/Qa,  ,  Q  are  varied  constants  and  a  is  the  angle 

of  incidence  with  0  a  =  30°. 

max 


NONLIFTING  VEHICLE  ENTRY  (a=0)  FROM  200  Km. 


INITIAL 

VELOCITY 

V 

0 

ENTRY 

ANGLE 

“*o 

LOADING 

(—2—) 

CDS 

ABSOLUTE 

ACCELERATION 

A 

max. 

TERMINAL 

CONDITION 

SEE 

FIGURES 

1.3911 

8° 

100  lb/ft2 

1 . 3g 

exit 

5,  11 

8-1/4 

100 

2.8 

exit 

11 

8-1/2 

100 

7.6}  7.0 

impact 

11 

9 

100 

18.1 

impact 

5,  11 

1.3911 

8 

25 

3.4}  8.6 

impact 

4,  12 

8-1/4 

25 

8.0 

impact 

4,  12 

8-1/2 

25 

13.6 

impact 

4,  12 

1.3088 

7-1/2 

100 

1.3 

exit 

15 

- 

8 

100 

8.9 

impact 

5,  15 

9 

100 

23.6 

impact 

5,  15 

1.3088 

7 

25 

1.02 

exit 

16 

7-1/2 

25 

3.3}  7.7 

impact 

16 

8 

25 

13.5 

impact 

16 

TABLE  I  "S"  Trajectories 
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LIFTING  VEHICLE  ENTRY  (V  =  1.3911)  FROM  200  km. 


INCI- 

)ENCE 

ANGLE 

a 

ENTRY 

ANGLE 

“Y0 

LOADING 

<!> 

DRAG 

at 

a=0 

C°o 

LIFT 

SLOPE 

9Cl 

0a 

LIFT 

DRAG 

& 

u 

max 

RELATIVE 

ACCELER¬ 

ATION 

T 

max 

TERMINAL 

CONDI¬ 

TION 

SEE 

FIGURES 

0 

8° 

25  lb/ft2 

0.6 

0.01 

1 

!  2.  Og 

exit 

6,8,13 

9  ' 

25 

0.6 

0.01 

1 

20.3 

impact 

6,7,13 

12 

25 

0.6 

0.01 

1 

53.8 

6,13 

8  * 

50 

0.4 

0.02 

0.5 

0.7 

exit 

14 

9 

50 

0.4 

0.02 

ra 

17.1 

impact 

7,14 

9 

50 

0.6 

0.01 

18.3 

impact 

7,14 

9 

25 

0.4 

0.02 

mm 

19. 1 

impact 

7,14 

5° 

8 

25 

0.6 

0.01 

i 

1.8 

exit 

8,17 

9 

25 

0.6 

0.01 

i 

15.0 

impact 

17 

10° 

8 

25 

0.6 

0.01 

i 

1.6 

exit 

8,17 

9 

25 

0.6 

0.01 

i 

1 2. 2 ; 5. 1 

impact 

17 

30° 

8 

25 

0.6 

0.01 

i 

1.6 

exit 

8,9,18 

9 

25 

0.6 

0.01 

i 

8.7 

exit 

9,18 

25 

0.6 

0.01 

i 

1 7 . 0  ;8 . 6  | 
3. 4  ;  1 . 2 

impact 

9,10,18 

25 

0.6 

0.01 

i 

31.0J1.8* 

1.1 

impact 

9,18 

TABLE  II  "B"  Trajectories 


A  complete  set  of  peak  temperatures  is  given  by' Tables 
IV  through  XII.  These  formed  the  basis  of  the  contours  of  Figs. 
4  through  18. 
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8 .  Temperature-Time  Histories 
Entry  without  lifti 

Figure  4.  The  effect  of  entry  angle  on  trajectory  (R)  , 
absolute  wall  temperature  (T)  ,  absolute  acceleration. 

VQ  =  1.3911,  CD?/m  =  0.04,  a  =  0,  r'  =  4,  D  =  0.01. 

Figure  5.  The  effect  of  initial  velocity  on  trajectory 
(R),  absolute  wall  temperature  (T)  ,  absolute  acceleration. 

-Y  =  8°  and  9°,  CpS/m  =  0.01,  a  =  0,  R  =  6,  D  =  0.04. 

Figure  6.  The  alleviating  effect  of  steep  entry  on 
peak  wall  temperature.  VQ  =  1.3911,  C^s/m  =  0.024,  a  =  0, 
r'  =  4,  D  =  0.03. 

Figure  7.  The  effect  of  aerodynamic  loading,  C^S/m. 

Vq  =  1.3911,  -y  =  9°,  a  =  0,  r'  =  6,  D  =  0.05. 


Entry  with  lift* 

Figure  8.  The  effect  of  constant  lift  on  wall  temperature. 

VQ  =  1.3911,  -y  =  8°,  m/S  =  25  lb/ft2,  CD  =  0.6,  0CL/aa  =  0.01, 

(C  /C  )  =  0.5,  r'  =  6,  D  =  0.01.  0 

max 

Figure  9.  The  effect  of  entry  angle  at  maximum  lift. 

Vq  -  1.3911,  a  =  30°,  m/S  =  25  lb/ft2,  CD  =  0.6,  aCL/a<x  =  0.01, 
(cl/cd)  =  0-5,  R'  =  6,  D  »  0.03.  ° 


max 


Figure  10.  Entry  at  maximum  lift}  the  effect  on  tempera- 

I  Q 

ture  of  nose  geometry  (R  and  D) .  VQ  =  1.3911,  a  =  30  , 

-Y_  =  10°,  m/S  =  25  lb/ft2,  Cn  =  0.6,  aC./aa  =  0.01, 

(CL/CD)  =  0.5.  0 


max 
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Influence  of  Nose  Geometry  on  Peak  Wall 


>ratur< 


Entry  without  lift! 

Figure  11.  The  effect  of  entry  angle.  VQ  =  1.3911, 
ci  —  0,  C^S/m  —  0.010. 

Figure  12.  The  effect  of  entry  angle.  VQ  =  1.3911, 
ci  —  0 ,  Cq  S  /m  —  0.040. 

Figure  13.  The  effect  of  steep  entry.  VQ  =  1.3911, 
ci  —  0,  CpS/m  =  0.024. 

Figure  14.  The  effect  of  aerodynamic  loading  C^S/m. 

V  =  1.3911,  a  =  0. 
o  7 

Figure  15.  The  effect  of  entry  angle.  VQ  =  1.3088, 
ci  —  0,  CpS/m  —  0.010. 

Figure  16.  The  effect  of  entry  angle.  VQ  =  1.3088, 
q  —  0,  CpS/m  —  0.040. 


Entry  with  lifti 

Figure  17.  Low  and  moderate  lift  entry.  VQ  =  1.3911, 

a  =  5°  and  10°,  -y  =  8°  and  9°. 

7  o 

Figure  18.  Maximum  lift  entry.  Vo  =  1.3911,  a  =  30°, 
-r0  =  8,  9,  10,  12°. 


10.  Conclusions 

a.  Temperature-time  histories,  and  particularly  peak  temperature, 
are  essentially  independent  of  the  entry  equilibrium  wall 
temperature  for  the  practical  range  of  entry  conditions. 

b.  For  multiple-pass  trajectories  the  vehicle  reaches  equilibrium 
temperature  during  the  ex-atmospheric  phase  prior  to  reentry 
unless  the  ex-atmospheric  trajectory  is  closely  circular  near 
the  sensible  limit  of  the  atmosphere  (100-120  km.  altitude). 

c.  Whereas  with  increasing  entry  angle  the  maximum  deceleration 
increases  rapidly,  the  peak  wall  temperature  is  fairly  in¬ 
sensitive  to  entry  angle,  and  moreover,  the  peak  value  rises 
to  a  maximum  and  then  decreases. 
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d.  Lower  initial  entry  velocities  at  the  same  entry  angle 
may  or  may  not  lower  the  peak  wall  temperature  and 
maximum  deceleration. 

e.  The  effect  of  decreasing  the  aerodynamic  loading  m/c^S 
lowers  the  peak  wall  temperature. 

f.  The  effect  of  entry  with  constant  lift  is  to  decrease 
slightly  the  peak  temperature,  monotoni cal 1 y  with  in¬ 
creasing  incidence  angle.  Peak  deceleration,  however, 
can  be  drastically  reduced  at  high  incidence. 

g.  "Skip"  trajectories  within  the  effective  atmosphere  can 
occur  with  or  without  lift.  These  are  characterized  by 
twin  deceleration  and  temperature  peaks  (for  a  nonlifting 
entry)  or  by  multiple  deceleration  peaks  for  a  lifting 
entry.  In  the  latter  case  there  may  be  fewer  temperature 
peaks  than  acceleration  peaks.  In  any  case  subsequent 
peak  values  of  temperature  decrease  in  magnitude,  although 
intervening  temperatures  usually  remain  moderately  high. 

Tne  first  deceleration  peak  may  or  may  not  be  critical; 
for  nonlifting  entry  the  second  peak  is  often  much  larger; 
for  lifting  entry  subsequent  deceleration  peaks  are  much 
smaller  than  the  first  for  the  cases  presented  here  but 
this  is  by  no  means  always  so. 

h.  Peak  temperatures  are  relatively  insensitive  to  the  vehicle 
nose  radius  by  comparison  with  the  effect  of  skin  thickness. 
The  higher  peaks  occur  of  course  with  the  smaller  thicknesses; 
but  the  latter  lead  to  lower  exit  or  impact  wall  temperatures 
when  heat  loss  due  to  reradiation  occurs. 

i.  The  curves  of  Figures  11-18  suggest  that  for  manned  or  un¬ 
manned  entry  to  impact  with  acceptable  decelerations  a 

I 

practical  combination  of  nose  radius  (R  )  and  skin  thickness 
(D)  exists.  For  example,  for  manned  entry  with  peak  deceler¬ 
ation  close  to  a  15g  limit  and  peak  wall  temperature  to 
0  * 

about  5000  R  reasonable  R  -D  combinations  can  be  found  for 
entry  under  the  conditions  given  by  Figs.  12-14,  16-18. 
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Of  these  perhaps  the  most  desirable  are  Figs.  16  and  18. 

The  former  is  a  non-lifting  ("S")  entry  from  VQ  =  1.3088 

with  an  8°  entry  angle  and  C^S/m  =  0.040  while  the  latter 

is  a  "B"  entry  with  maximum  lift  from  VQ  =  1.3911  with  a 

10°  entry  angle  and  CD  S/m  =  0.036. 

o 
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APPENDIX 


Ex-Atmospheric  Temperature  Response  After  Exit 


Figures  1,  2  and  3  enable  the  determination  of  the  time- 
temperature  or  the  position-temperature  relationship  of  the 
vehicle  during  ex-atmospheric  motion  between  exit  from  and  re¬ 
entry  into  the  planetary  atmosphere. 

At  exit  the  following  quantities  are  requiredt 
e  eccentricity 

P  nondimensionalized  semi-latus  rectum,  p/r 

O 

u  true  anomaly 

x 

T  temperature  at  exit 

R  exit  altitude  (nondimensionalized) 

These  quantities  fix  the  ex-atmospheric  flight  time  and  so, 

for  a  given  vehicle,  the  reentry  temperature.  To  find  the 
latter  the  following  procedure  can  be  used* 

(i)  Compute  the  nondimensional  semi-major  axis  A  =  p/( 1 -e  ). 
Use  Fig.  1  to  find  t  the  period  of  the  ellipse. 

(ii)  From  Fig.  2  find  the  nondimensional  time  of  flight 
At  /t  using  the  exit  value  of  u  together  with  e. 

(iii)  Then  At  =  (At  /t)t.  „ 

(iv)  In  Fig.  3  locate  the  time  +  1000)  corresponding 

x  u  " 

to  T  and  scale  it  for  wall  thickness  D  and  so  find  t  . 

I 

(v)  To  this  time  add  At  .  Normalize  the  resultant  time  to 

M  . 

D  =  0.02,  i.e.,  recalculate  (0.02t  /D  +  1000)  and  read 
from  Fig.  ^  the  reentry  temperature.  See  also  Table  III. 


I 


T[°R] 


t "  ( T )  (sec.) 


t " ( T )  (sec.) 
0 


711.13 

+  oo 

1700 

711.2 

30,034 

1710 

711.5 

24,281 

1725 

712 

21,328 

1750 

720 

13,318 

1800 

800 

5604.6 

1900 

900 

3331.4 

2000 

1000 

2196.4 

2130 

1100 

1497.0 

2200 

1200 

1024. 5 

2300 

1300 

688.39 

2400 

1400 

440.95 

2500 

1500 

254.40 

3000 

1600 

111.29 

3550 

1650 

52.223 

4000 

1675 

25.316 

4500 

1690 

9.9375 

5000 

1700 

0 

71,150 

D  =  0.02  feet 


-  9.6250 

-  23.781 

-  46.422 

-  88.203 
-159.88 
-218.86 
-281.0 
-309.08 
-343.98 
-373.80 
-399.45 
-486.  i9 
-537.47 
-562.80 
-581.48 
-594.25 
-640.34 


TABLE  III 
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V  =  1.3911 


m 

CD  S/m 
0 

D  ft. 

(dl) 

'  dt' 

max 

Approx. 
At^  sec 

T  °R 

max 

■ 

g 

3max 

TERMINAL 

CONDITION 

8 

0.01 

9 

0.119 

60 

6000 

-10 

1.3 

EXIT 

9 

f 

0.0889 

70 

5500 

-40 

9 

0.01 

0.105 

55 

5490 

-10 

0.0735 

70 

4840 

-50 

9 

0.04 

0.0548 

75 

4090 

-65 

8.25 

0.01 

9 

gg 

0.142 

55 

m 

5 

2.8 

EXIT 

9 

0.114 

60 

H9 

-25 

9 

0.01 

0.128 

60 

6130 

-5 

0.03 

0.0964 

70 

5640 

-35 

9 

0.04 

0.0744 

75 

4940 

-50 

8.5 

0.01 

9 

0.01 

0.165 

55 

7200 

m 

IMPACT 

0.142 

45 

6780 

■ 

9 

0.151 

55 

6540 

B 

0.03 

0.121 

55 

6040 

9 

0.04 

0.0959 

60 

5240 

-10 

B 

9 

0.01 

9 

0.01 

0.214 

30 

7630 

18.1 

IMPACT 

0.196 

30 

7000 

0.194 

30  ! 

6950 

5 

0.03 

0.164 

25 

6070 

5 

91 

0.04 

0.132 

30 

5060 

-5 

"S"  TRAJECTORIES 

At,  =  [Time  at  T  ]  -  [Time  at  (^r)  ] 

1  max  d  t 

max 


At0  =  [Time  at  g  -  Time  at  T  ] 

2  smax  max 


TABLE  IV 
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V  =  1.3911 


1 

o 

o 

o 

CO, 

3 

9 

(«) 
dt  max 

■ 

T 

ma 

1st 

PE 

°R 

X 

'  2nd 
IK 

Approx. 

sec 

1st 

PE 

1 

TERM¬ 

INAL 

CONDI¬ 

TION 

8 

0.04 

2 

1  0.01 

0.123 

5650 

438C 

-5 

3.4 

8.6 

IMPACT 

2 

0.101 

5420 

424C 

-15 

4 

0. 105 

60 

5170 

401C 

-5 

2 

0.03 

0.086 

70 

5080 

406C 

-35 

3 

0.02 

0.0920 

60 

5070 

400C 

-25 

4 

0.03 

0.0686 

80 

4460 

362C 

-45 

8.25 

0.04 

Bj 

0.01 

0.136 

40 

60' 

10 

■9 

B 

2 

IMPACT 

3 

0.02 

0.121 

40 

5720 

119 

■ 

B 

0.01 

0.123 

40 

■ 

80 

■ 

B 

0.02 

0.110 

50 

3 

60 

■ 

B 

0.03 

0. 102 

50 

5300 

50 

■ 

9 

0.03 

0.0846 

50 

4620  ! 

50 

■ 

8.5  : 

0.04 

2 

0.01 

0.160 

40 

6230 

10 

13. 

6 

IMPACT 

2 

0.02 

0.142 

40 

5810 

-5 

4 

0.01 

X).  142 

30 

5680  ■ 

10 

3 

0.02 

0.127 

40  ! 

5420  ‘ 

-5 

2 

0.03 

0.122 

40 

5280 

-15 

4 

0.03 

0.101 

40 

4560  , 

-15  . 

TABLE  V 
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V  =  1.3088 


-Y0deg. 

o 

o 

o 

cn, 

3 

• 

R 

ft. 

D  ft. 

(dT) 

'  d  v 

max 

Approx. 
Atj  sec 

T  °R 

max 

Approx. 
At2  sec 

q 

3max 

TERMINAL 

CONDITION 

7.5 

0.01 

2 

0.01 

0.0995 

60 

5770 

-5 

1.3 

EXIT 

2 

0.03 

0.0759 

85 

5330 

-40 

4 

0.01 

0.0889 

70 

5280 

-10 

4 

0.03 

0.0627 

75 

4710 

-50 

6 

0.04 

0.0470 

90 

4000 

-70 

8 

0.01 

2 

0.01 

0.144 

55 

6920 

35 

8.9 

IMPACT 

2 

0.03 

0.126 

45 

6460 

20 

4 

0.01 

0.132 

45 

6340 

35 

4 

0.03 

0.107 

45 

5730 

15 

6 

0.04 

0.0842 

60 

4920 

-5 

9 

0.01 

2 

0.01 

0.229 

25 

7500 

20 

23.6 

IMPACT 

4 

0.01 

0.216 

25 

6820 

20 

2 

0.03 

0.202 

30 

6430 

5 

4 

0.03 

0.175 

35  | 

5590 

5 

6 

0.04 

0.137 

30 

4520 

-5 

TABLE  VI 
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V  =  1.3088 


-YQdeg. 

CD  S/n 
0 

t 

R 

ft. 

D  ft. 

(dl) 

vdt; 

max 

Approx. 
At^  sec 

T  °R 

max 

Approx. 
At2  sec 

g 

3max 

TERMINAL 

CONDITION 

7 

0.04 

2 

0.01 

0.0610 

70 

4180 

-25 

1.0 

EXIT 

3 

0.01 

0.0550 

50 

3920 

-25 

4 

0.01 

0.0520 

60 

3740 

-35 

2 

0.02 

0.0466 

80 

3710 

-55 

5 

0.01 

0.0494 

60 

3610 

-35 

3 

0.02 

0.0407 

70 

3410 

-55 

2 

0.03 

0.0365 

90 

3260 

-75 

2 

0.04 

0.0297 

90 

2900 

-75 

4 

0.03 

0.0282 

90 

2760 

-75 

7.5 

0.04 

2 

0.01 

0.0808 

80 

5430 

-5 

3.3* 

IMPACT 

2 

0.02 

0.0868 

60 

5180 

-25  | 

4 

0.01 

0.0792 

40 

4950 

-5 

2 

0.03 

0.0729 

80 

4850 

-45 

3 

0.02 

0.0789 

60 

4850 

-25 

4 

0.03 

0.0584 

90 

4250 

-65 

8  ! 

0.04 

2 

0.01 

0.139 

30 

5940 

20 

13.5 

IMPACT 

2 

0.02 

0.124 

30 

5480 

10 

4 

0.01 

0.127 

40 

5370 

10 

3 

0.02 

0.114 

40 

5080 

-5 

2 

0.03 

0.106 

40 

4960 

-5 

4 

0.03 

0.0856 

40 

4240 

-15 

# 

second  g 

3max 


7.7 


TABLE  VII 
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V  =  1 

3911 

p 

11 

0 

C^S/m 

=  0.024 

D  ft. 

(dT) 

vdt' 

max 

Approx . 
Atj  sec 

T  °R 

max 

Approx. 
Atj  sec 

q 

3max 

TERMINAL 

CONDITION 

8 

D 

0.01 

0.118 

55 

5840 

0 

2.0 

EXIT 

0.0964 

;  50 

5050 

-10 

B 

0.0653 

80 

4660 

-55 

0.03 

0.0702 

65 

4650 

-40 

B 

0.05 

0.0446 

90 

3620 

-75 

9 

B 

0.01 

0.210 

30 

6930 

mm 

20.3 

IMPACT 

0.180 

25 

5920 

B 

0.141 

30 

4990 

0.05 

0.132 

30 

4790 

-10 

B 

0.05 

0.0914 

35 

3580 

-15 

12 

B 

0.01 

0.450 

1 

7370 

0 

53.8 

IMPACT 

0.376 

6060 

0 

B 

0.316 

10 

5260 

0.256 

15 

4380 

0.05 

0.232 

15 

4060 

-10 

B 

0.05 

0.157 

15 

2960  ! 

-10 

"B"  Trajectories 


H  T 

At.  =  [Time  at  T  ]  -  [Time  at  (-77)  ] 

1  u  max  dt 

max 

m2  =  CTlme  at  gmax]  -  [Tima  at  T„,ax] 


TABLE  VIII 
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V  =  1 . 3911  a  =  0 


HjEffl 

CD  S/m 
0 

D  ft. 

(dl) 

^dt' 

max 

Approx . 
At2  sec 

T  °R 

max 

Approx. 

sec 

g 

3max 

8 

0.008 

12 

0.01 

0.119 

60 

ESI 

-10 

0.7 

EXIT 

6 

0.01 

0.0975 

55 

-15 

2 

0.05 

0.0687 

75 

-60 

4 

0.03 

0.0734 

70 

EH1' 

-50 

6 

0.05 

0.0475 

80 

3740 

-70 

9 

0.008 

9 

0.01 

0.215 

35 

7870 

15 

17.1 

IMPACT 

9 

0.01 

0.190 

30 

6830 

12 

9 

0.03 

0.172 

30 

j  6340 

0 

9 

0.05 

0.163 

35 

6320 

-5 

9 

0.05 

0.121 

35 

4890 

-10 

9 

0.012 

2 

mm 

0.215 

30 

7520 

15 

18.3 

IMPACT 

6 

0.189 

25 

6500 

10 

4 

yufl 

0.162 

30 

5840 

0 

2 

0.05 

0.153 

35 

5740 

-5 

5 

0.04 

0.132 

30 

5040 

-5 

6 

0.05 

0.110 

35 

4380 

-15 

9 

0.016 

2 

0.01 

0.212 

>20 

>7170 

<20 

19.2 

IMPACT 

4 

0.01 

0.197 

30 

6640 

10  ! 

3 

0.02 

0.189 

25 

6420 

5 

2 

0.03 

0.182 

25 

6340 

0 

6 

0  01 

0.184 

25 

6270 

10 

3.5 

0.03 

0.159 

30 

5650 

-5 

4 

0,03 

0.153 

>25 

>5380 

<5 

2 

0.05 

0.145 

30 

5340 

-5 

6 

0.03 

0.137 

30 

5010 

-5 

6 

0.05 

0.102 

35 

4040 

-15 

TABLE  IX 
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V  =  1.3911  a  =  30 


M 

D  S/m 
0 

m 

1  t  t  • 

<3T> 

max 

Approx. 
At^  sec 

T  °R 

max 

Approx. 
At2  sec 

^max 

TERMINAL 

CONDITION 

8 

0.024 

2 

0.01 

0.114 

40 

EH 

-10 

1.6 

EXIT 

4 

0.01 

0.103 

45 

-15 

3 

0.02 

0.0872 

55 

ESI, 

-30 

6 

0.01 

0.0938 

50 

4670 

-15 

2 

0.03 

0.0807 

55 

4630 

-40 

3.  5 

0.03 

0.0678 

65 

4090 

-50 

4 

0.03 

0.0648 

65 

3980 

-50 

2 

0.05 

0.0590 

80 

3850 

-65 

6 

0.03 

0.564 

65 

3620 

-50 

6 

0.05 

0.0393 

90 

2890 

-80 

9 

0.024 

2 

1 

0.205  i 

EH 

0 

8.7 

EXIT 

4 

0.152 

-15 

2 

0.129 

40 

6 

0.02 

0.136 

30 

H 

2 

0.05 

0.112 

50 

ESS 

-40  ; 

6 

0.03 

0.108  | 

45 

4130  ! 

-35 

4.5 

0.05 

0.0844 

60 

3570 

-40 

TABLE  XI 
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V  =  1.3911  a  =  30 


FIRST 

PEAK 

SECOND  PEAK 

-*0 

deg 

CD  S/m 
0 

f 

R 

ft. 

D  ft 

(dl) 

max 

.  Appx 
'  At^ 

sec 

T  R° 

ma  x 

Appx 

At2 

sec 

<3T> 

max 

Appx 

Atl 

sec 

T 

max 

Appx 

*t2 

sec 

TERMIN¬ 

AL 

CONDI¬ 

TION 

10 

0.024 

;  2 

0.01 

0.287 

20 

6740 

0 

0.0818 

20 

4230 

-5 

IMPACT 

6 

0.01 

0.238 

15 

5640 

-5 

0.0628 

25 

3490 

-10 

;  4 

0.02 

0. 206 

25 

5160 

-15 

0.0434 

25 

3290 

-15 

I  6 

0.02 

0.183 

25 

4720 

-15 

0.0373 

25 

3070 

-15 

!  2 

0.04 

0.172 

35 

4700 

-25 

0.0308 

35 

3350 

-25 

4 

0.03 

0.164 

35 

4460 

-25 

0.0304 

35 

3100 

-25 

2 

0.05 

0.149 

45 

4250 

-35 

0.0248 

35 

3280 

-25 

6 

0.03 

0.144 

35 

4040 

-25 

0.0257 

35 

2920 

-25 

4.5 

0.05 

0.112 

45 

3440 

-35 

0.0178 

35 

2900 

-25 

6 

0.05 

0.100 

45 

3180 

-35 

0.0158 

35 

3100 

-25 

•i.x  -  n-° 

g2 

ymax 

=  8.6 

g3  =  3.4 

ymax 
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FIG.  I  MAJOR  AXIS  &  PERIOD  OF  KEPLERIAN  ELLIPSE 


NONDIMENSIONAL  EX -ATMOSPHERIC  TIME  (At'/r') 


FIG.  2  EX-ATMOSPHERIC  FLIGHT  TIME  IN  TERMS  OF  EXIT  CONDITIONS 
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FIG.  3  REENTRY  TEMPERATURE  AND  EX -ATMOSPHERIC  FLIGHT  TIME 


TIME  (SECONDS) 

FIG.  4  EFFECT  OF  ENTRY  ANGLE  ON  WALL 
TEMPERATURE  AND  DECELERATION 
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FIG.  5  EFFECT  OF  ENTRY  VELOCITY  ON  WALL 
TEMPERATURE  AND  DECELERATION 


I 


FIG  6  ALLEVIATING  EFFECT  OF  STEEP  ENTRY 
ON  WALL  TEMPERATURE 


1000  1500 
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FIG.  9  EFFECT  OF  ENTRY  ANCLE  AT  MAXIMUM  LIFT 
ON  DECELERATION  AND  TEMPERATURE 
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(a)  THE  EFFECT  ON  WALL  TEMPERATURE  OF  NOSE 

SKIN  THICKNESS 


TIME  (SECONDS) 

(b)  THE  EFFECT  ON  WALL  TEMPERATURE  OF  NOSE  RADIUS 

FIG.  10  ENTRY  AT  MAXIMUM  LIFT: 

WALL  TEMPERATURE  -  TIME  HISTORIES 
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FIG.  II 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
EFFECT  OF  ENTRY  ANGLE  AT  ZERO  LIFT 


NOSE  RADIUS  R’  (FEET) 

V0=l.39ll  CpS/m  =0.040  FT^LB.  or  -  0 


FIG.  12 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
EFFECT  OF  ENTRY  ANGLE  AT  ZERO  LIFT 
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FIG.  13 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
EFFECT  OF  STEEP  ENTRY  AT  ZERO  LIFT 
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FIG.  14 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
EFFECT  OF  m/CoS  AT  ZERO  LIFT 
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V0- 1.3088  CoS/m  -  0.010  FT*/|_B.  a-0 


FIG.  15 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
EFFECT  OF  ENTRY  ANGLE  AT  ZERO  LIFT 
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FIG.  16 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
EFFECT  OF  ENTRY  ANGLE  AT  ZERO  LIFT 
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FIG.  17 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE 

LOW  LIFT  ENTRY 
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FIG.  18 

NOSE  GEOMETRY  AND  PEAK  WALL  TEMPERATURE: 
MAXIMUM  LIFT  ENTRY 


